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Abstract 

Recently obtained NNLO exact corrections for Higgs and Pseudoscalar Higgs 
boson production in hadron colliders are compared with approximate ones. As shown 
before, it is found that there is a range of a proper variable where these corrections 
differ little. 



*e-niail: apcont@physics.nicgill.ca, acontog@cc.uoa.gr 
^e-mail: ppapachr@cc.uoa.gr 



Some time ago it was argued that for processes involving structure functions and/or 
fragmentation functions, over a range of a proper kinematic variable w, there is a part that 
dominates the next- to leading order (NLO) correction and that this part contains the 
distributions 5(1 — w) and [ln"(l — w)/l — w)]+ = 0, 1 p. Subsequently this argument 
was extended to the then existing next- to- next- to leading order (NNLO) calculations, 
namely Drell- Yan (D-Y) production of lepton pairs {q + q 7*) and deep inelastic 
structure (DIS) functions {q + l* ^ q) [2 *. 

In the meantime two more processes have been calculated in NNLO: Higgs boson pro- 
duction in hadron- hadron collisions [g + g H) |3j and neutral pseudoscalar Higgs boson 
production in hadron- hadron collisions {g + g ^ A) Clearly, it would be important 
to see whether the procedures developed in |2] apply also to Higgs and pseudoscalar Higgs 
boson production as well. 

In the calculation of Higgs boson production in hadron- hadron collisions to leading 
order (LO) |3] and to NLO "Bj no approximations of the Higgs- two gluon vertex are 
necessary''". This vertex is dominated by the top quark, which is known to have a mass mt 
much greater than that of the other quarks. However, the NNLO calculation was possible 
only in the limit of the Higgs mass mu 

iJiH « 2mt. (1) 

In this limit the top-quark loops are replaced by point-like vertices and the corresponding 
effective Lagrangian [HI [3 El is known to provide a satisfactory description of the cross 
section for a Higgs boson at NLO [H]. 

In the calculation of the pseudoscalar Higgs boson production the situation is more 
complicated. The Higgs boson sector of the Minimal Supersymmetric Standard Model 
consists of two complex Higgs doublets. Thus, apart from the mass of the neutral pseu- 
doscalar Higgs boson m^, the ratio of the vacum expectation values of the two Higgs 
doublets 1^1/^2 = tan/3 also enters. The calculation of |4j is valid for small and moderate 
values of tan/3; only then the gg ^ Ais dominated by a top-quark loop. Then, for 

tra « '2mt (2) 

the interaction of the pseudoscalar Higgs boson can be described by an effective Lagrangian 

m 

In our approach j2] the proper variable is proportional to r = rri\j/S (or r = m\/S), 
where is the total cm. energy of the initial hadrons. Our approach requires inclusion 
of the region r large (inclusion of r near 1). Since experiment excludes values of ttih (or 
t^a) < 100 GeV, we have to consider \/S well exceeding this value. On the other hand, 
inclusion of \^ > 2 TeV would require tuh (or m^) well exceeding 1 TeV, which would 

*This paper, in preliminary form, was presented in "QCD 2000", Montpellier, France, July 2000, and 
published in Nucl. Phys. B, (Proceed. Suppl.) 96 (2001) 94 fhep-ph/0109070). 

^To NLO this is due to an accidental cancellation of the dependence on the top quark mass between 
real and virtual corrections. See the last paper of Ref. [6]. 
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render questionable the field- theoretic approach. We then have considered a nominal en- 
ergy of = 520 GeV. Clearly, for tuh (or m^) > 200 GeV the inequality (1) (or (2) ) 
is violated and the whole results of J3| and p[J, which we use, should be considered as just 
providing a mathematical model, where the approach of can be tested. 

We begin with pp ^ H + X (or pp — > if + X) mediated via the subprocess gg H* 
and we consider the cross- section 

1 

(Th^+h2->H+x{mjj,S) = J dXidX2fg/p{xi)fg/p{x2)(Tgg^H{mH,XiX2S) (3) 



where hi, h2 denote p,p (or p,p) and fg/p{x) is the standard distribution of gluons inside 
the p (or p). Using dimensional analysis we write the partonic cross- section in terms of 
the dimensionless variable 

(4) 



my T 



X1X2S X1X2 



and after factoring the coUinear singularities (usually in the MS scheme) we end up with 
the following expression [S] 

ah,+h2^H+xiT, S) = rfg/p ® fg/p ® {agg{z)/z) (r) (5) 

where ® denotes the standard convolution defined as 

1 

[/l ® /2] (r) = / dXidX2fl{Xi)f2{x2)S{T - X1X2). (6) 



The partonic cross- section (ygg{z) is given by the following perturbation expansion 

(7) 
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where the functions rjf^, = 0, 1, 2, are given in Eqs. (44), (45), (47), (48) and (49) of 0, 
and 

TT /as 
° 576t;^ V vr 
with V ~ 246 GeV the Higgs vacum expectation value. 



Subsequently we proceed as in |2]^ We write, for simplicity, (T/ji+/j2^h+x(t, 5*) = o-//(r, 5*) 
and denote by <J^h\t, S), k = 0, 1, 2, the 0{a^) part of cr//(r, S), by a^] the part of cr^'' 

*We remind that in our approach QIE] the various perturbation orders (LO, NLO, NNLO) should refer 
to the same subprocess. 

^Although known since long ago (see [1]), as in [2], we present also results for k = 1. 
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arising from distributions 6{1 — z) and [ln"(l — z)/l — 2;)]+ , here n = 0, 1,2,3 (virtual, 
collinear and soft gluons) and by a^^l the rest. We also define 



L^Pir,S) = ^-J^. (9) 



In the subsequent calculations we use Uf = 5 flavors and fix the renormalization and 
factorization scales at /i = M = rriH- For the gluon distributions we use the updated MS 
CTEQ5M1 set of 

Fig. 1, upper part, shows L^h\ A; = 1, 2, as functions of ^/T. For L^h\ while for relatively 
small i/r is significant, for ^/T > .63 is below 30%. As for L^^\ for ^/T > .43 is smaller 
than 20%. Moreover, both L^"* decrease fast as ^/T increases towards 1. 

As in P^, it is of interest to see the precentage of cr^^ of the total cross section de- 
termined up to 0(aj). Fig. 1, upper part, also shows the ratios (7^1,/ (o'^h^ + (^h^^ ^-^^ 
'^Hh/ {^H^ + '^H^ + '^H^\ The former is below 31% and the latter below 15% for all \^*. 
Again, both ratios decrease fast as y/r increase. 

Now we turn to the calculation of the pseudoscalar Higgs boson production and consider 
pp ^ A + X (or — > A + X) mediated via the subprocess gg — > A. As before, the partonic 
cross-sections crgg{z) have an expansion similar to (7) 



(10) 



where z is given by (4) with r = m\/S, (t>^gg{z) are given in Eqs. (8)-(ll) of ^ (together 
with the expressions of ri'^^g{z))^ and here 

^° = 256t;nan2/3lvj ' ^^^^ 

Now we write ia.T? j3ah^+h2->A+x = o"a(t, S) and, as before, denote by cr^'*(T, S) the 0(aj) 
part of crA(T, S), by a^j the part of cr^^ arising from distributions and by cr^^ the rest. We 
define 

Ly(r,5) = 4^ (12) 
crA{r,S) 

and fix the renormalization and factorization scales at /i = M = rriA- Again, for the gluon 
distributions we use the set CTEQ5M1 set of |10^. 

Fig. 1, lower part, shows L^a\ k = 1,2, as functions of ^/T. All the results are 
similar as for L^"*. Similar are also the results for the ratios cr^Ah/ ('^A^ + '^a'') 

*This can also be seen in the first paper of Ref. [11]. 
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We note the following*: suppose that in aj^] and a^fj, apart from the terms arising from 
the distributions 5(1 — z) and [ln"'(l — z)/l — 2;)]+ we include also the terms ln™'(l — z), 

m = 1,2, 3. Defining as a^^l and a^^l the rest, we find that the ratios u'^Hh/u^H CF^Ah/'^^A 
decrease significantly in magnitude over the entire range of a/t. Of course, the same holds 
for the ratios cr^)^/ [y'P + <7^h^ , o"!/!/ (c"!/'' + <^^h + ^h'*) ^^^1 the corresponding ratios with 
H replaced by A. 

Note also that R.ef . [TT] has obtained numerical results very similar to ^ and by ex- 
panding the phase-space integrals around the kinematic point z = t/xiX2 = 1, where 
T = m'jj/S or r = m\/S, and keeping a number of terms. Although the first paper of [11] 
was pubhshed before [3], we prefer the methods of [3] as they avoid expansions.''" 

Finally, in Fig. 2, upper part, we present the total cross-sections a'j^ + a'^j + (dashed 
line) and cr^^ + a^^^ + a^^^ (solid line). What is important is that as ^JT increases towards 
1 both cross-sections approach each other , and for r > 0.8 practically coincide. The same 
is observed in Fig. 2, lower part, which shows the quantities a^^ (dashed) and 

<yf + <^^ll + <yfs (solid). 

In conclusion, under the assumptions discussed at the beginning, we have shown that 
not only in D-Y production and DIS |2j, but also in Higgs and pseudoscalar Higgs boson 
production {gg —>■ H and gg — A) there is a part containing the distributions 5(1 — z) 
and [ln"'(l — z)/l — z)]+, here n = 0,1,2,3 (virtual, soft and coUinear part) that for 
A/r(= rriij/y/S or niA/y/S) not too small dominates the NLO and NNLO correction. This 
part is determined much easier than the NLO and in particular the NNLO correction. Of 
course, as it was stressed in [2], this part should not be restricted to too small a region 
near 1, for threshhold resummation ^12j becomes very important. 



NOTE ADDED 

After the completion of this article, the paper by V. Ravindran, J. Smith and W. van 
Neerven, hep-ph/0302335, appeared, confirming the results of [3], [4] and [11] by a differ- 
ent method. 
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Figure 1: Upper part: The ratios L^^ and a^^hl i^^^H + ^i/^) (dashed hnes) and the 
ratios L^^ and ct^hhI {^^^h ^^h ^^h^ (sohd hnes) versus = uih/VS- Lower 
part: The quantities and o-^aI/ {'^a^ + ^a"*) (dashed) and the quantities L^"* and 
^Ah/ (^a'' + ^A^ + (solid) versus ^/T = VHaI^- 
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Figure 2: Upper part: The cross-sections a\j +(7^^ +a\j (dashed Une) and cth +o'hs~'<~<^hs 
(sohd hne) versus a/t = mul^fS. Lower part: The quantities cr^'' + + cr^'' (dashed) 
and cr^^ + cr^^j + cr^^j (sohd) versus ^Jr = tua/VS. 
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